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We investigate the use of a non-degenerate parametric oscillator (NDPO) as a source for quantum 
lithography, for which the light can have high-flux and strong non-classical features. This builds on 
the proposal of Boto, et al. [A. N. Boto, et al., PRL 85, 2733 (2000)], for etching simple patterns on 
multi-photon absorbing materials with sub-Rayleigh resolution, using special two-mode entangled 
states of light. An NDPO has two outgoing modes differentiated by polarization or direction of 
propagation, but sharing the same optical frequency. We derive analytical expressions for the 
multi-photon absorption rates when the NDPO is operated below, near, and above its threshold. 
The resulting interference patterns are characterized by an effective wavelength half that for the 
illuminating modes. We compare our results with those for the case of a high-gain optical amplifier 
source, and discuss the relative merit of the NDPO. 

PACS numbers: 42.50.-p, 42.65.Yj, 42.50.St 
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I. INTRODUCTION 

The Rayleigh criterion states that diffraction limits the 
resolution of a traditional-optical lithographic system, 
and that the minimum feature size is determined by half 
the wavelength A of the illuminating beam. This limit 
arises in part from the fundamental photon statistics of 
laser light, according to which the constituent photons 
behave as if they are uncorrelated. In the year 2000, 
Boto, et al., set out a concrete proposal for exploiting 
non-classical states of light to achieve resolution beyond 
the classical limit (termed super-resolution) in a lithog- 
raphy experiment [l[ . Their proposal has attracted con- 
siderable research interest 

H H EL H IS 0, 1 • A key idea 
is to use a source of path-entangled states of light of the 
form \N0) + \0N), in the photon-number basis, popularly 
termed "NOON" states. In the scheme, the NOON states 
propagate through a simple interferometer, and then in- 
terfere at a N-photon-absorbing recording material. In 
principle, the procedure can etch a series of straight lines 
corresponding to an effective wavelength A/ (27V). Varia- 
tions of the method have been proposed for creating more 
general one- and two-dimensional interference patterns, 
by employing a family of entangled states However, 
considerable work remains to be done to demonstrate 
the feasibility of the method. A general program aims to 
do so, and requires investigations in turn of the source, 
imaging system, and multi-photon absorption process. 

For the case of two-photon quantum lithography, one 
could use parametric down-conversion in a medium ex- 
hibiting a x' 2 ) optical non-linearity as a source of photon 
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pairs. Due to the Hong-Ou-Mandel effect, a pair of indis- 
tinguishable photons simultaneously incident on a sym- 
metric beam splitter will yield a 2-photon N00N-st&te 
at the output ports [l(|. Subsequent quantum interfer- 
ence of the two modes on a 2-photon absorbing mate- 
rial will then create a fringe pattern with half the fringe 
spacing achieved using a classical light source. However, 
the output of common parametric down-conversion ex- 
periments is not expected to be sufficiently bright to 
be useful for typical 2-photon recording materials. In 
Ref. 0, Agarwal, et al., considered a strongly-pumped 
high-gain optical parametric amplifier. By assuming a 
single-mode operation, it can be shown that the state 
generated by an unseeded optical parametric amplifier 
(OPA) is the two-mode squeezed vacuum state, which has 
the form X^^Lo tanh n (G)|n)|n) in the basis of photon- 
number states. Here, G is a gain parameter, which de- 
pends on the interaction volume in the crystal, the am- 
plitude of the electric field of the pump beam, and the 
strength of the second-order susceptibility \^ ■ To quan- 
tify the contrast of the interference pattern at the output, 
the visibility is defined as the difference of the maximum 
and minimum absorption rates, divided by sum of these 
rates. The visibility varies between and 1. Explicit 
calculation reveals that the visibility falls from 1 to an 
asymptotic value of 0.2, as the gain parameter G of the 
OPA is increased from — ► oo. This result was general- 
ized in Ref. Q to the case of higher-order multi-photon 
absorbing materials. Essentially the same behavior is 
seen in these cases, with a halving of the fringe spac- 
ing, and the visibility falling from 1 to an asymptotic 
value, which is greater the higher the order of the absorp- 
tion process. These predicted features have been demon- 
strated experimentally, using coincidence measurements 
at photodetectors to simulate the recording medium [8| . 

In this paper, we focus on using a NDPO (non- 



2 



degenerate optical parametric oscillator) as a practical 
high-flux source of non-classical light for quantum lithog- 
raphy. In this case, the process of parametric ampli- 
fication occurs in an optical cavity, in resonant signal 
and idler modes. The signal and idler modes share the 
same optical frequency, but have orthogonal polariza- 
tion. The cavity modes are coupled to external prop- 
agating modes through a transmissive end-mirror. Pho- 
tons are created in pairs in the intracavity modes, but 
evolve independently out of the cavity on a timescale 
characterized by the cavity lifetime. The twin beams 
at the output are highly correlated. Unlike the OPA, 
the NDPO has a well-defined threshold for oscillation, 
and the below-, near- and above-threshold regimes re- 
quire different mathematical treatments. The theory 
for the NDPO is complicated by the need to account 
for mode losses, and the coupling of cavity and exter- 
nal modes. To analyze the quantum lithography proce- 
dure, we will apply a methodology developed and ap- 
plied in a series of papers which investigate the proper- 
ties of the output fields of the optical parametric oscilla- 
tors Ell EH 0, EH . By using the positive-P representa- 
tion [151 ] , the Markovian master equation describing the 
dynamics of the standard model of the NDPO can be 
translated into a multi-variate Fokker-Planck equation. 
The Fokker-Planck equation can be treated using tech- 
niques from nonequilibrium classical statistical mechan- 
ics [H, E3], allowing steady-state expectation values of 
various observables to be evaluated. 

In Sec. UH we describe the setup for a quantum lithog- 
raphy experiment, and present a quantum-classical cor- 
respondence based on the positive-P representation for 
the master equation for the density matrix. A solution 
of the resulting Fokker-Planck and Langevin equations 
is presented following Refs. [H, EH- ^ n Sec. EE wc 
propagate the c-number variables for the NDPO's out- 
put modes through a simple interferometer, and deter- 
mine the multi-photon absorption of arbitrary order. We 
then look in detail at the features of the interference pat- 
terns, focusing on the below-threshold case in Sec. IIII Al 
and the near- and above-threshold case in Sec. IIII Bl In 
Sec. IIVI we compare our results with those for the OPA, 
and comment on their similarities and differences. 



II. METHOD 

A. The setup 

The setup for a quantum lithography experiment is il- 
lustrated in Fig. [TJ We start at the NDPO source. The 
modes corresponding to the signal and idler have fre- 
quency u>, and we denote the corresponding creation op- 
erators as a\ and a\. The signal and idler modes are 
assumed to experience linear losses characterized by the 
decay rate 7, which arises from transmission, absorption, 
and scattering losses. We assume the pump field to be 
a classical field of frequency 2lo and of normalized am- 




FIG. 1: A strongly-pumped NDPO serves as the source, with 
signal, idler and pump modes interacting via a — non- 
linearity inside an optical cavity. Signal and idler modes are 
frequency degenerate and have different polarizations. Exter- 
nal signal and idler modes are separated into different spatial 
paths by a polarizing beam splitter. The two modes are then 
combined at a 50:50 beam splitter, and are recombined at a 
multi-photon-absorbing recording material. 



plitude denoted by e (assumed to be given by a positive 
value). The cavity pump mode is subject to linear loss 
with decay rate 73, and we assume that 73 3> 7. The cre- 
ation operator corresponding to the cavity pump mode 
is S3. The Hamiltonian describing this system in the in- 
teraction picture can be written as, 



Hi = ihK [Ji\ala3-aia2alj + ih-f 3 e (al-a 3 J + Pq oss , 

where k is the mode-coupling constant determined by the 
strength of the optical nonlinearity, and H\ oss describes 
the interaction of the principal modes with all the reser- 
voir modes. Light emitted by the cavity into the external 
signal and idler modes, assumed to have orthogonal po- 
larizations, is separated into two spatial paths by a polar- 
izing beam splitter. These modes are then combined at 
a symmetric beam splitter, and finally recombined at the 
multi-photon absorbing material using mirrors, such that 
the beams are counterpropagating at grazing incidence 
over the area of interest. Two classical beams incident on 
single-photon absorbing material would generate a fringe 
pattern of the form cx [l + cos(2A;x)], where k denotes 
the optical wave number for the signal and idler modes, 
and x denotes translation across the substrate. We shall 
denote the optical phase difference 2kx by <j>. When the 
recording medium is a p-photon absorber, and the illumi- 
nating beams have statistics which need not be classical, 

the absorption rate is given by er^tr (j)E + (xf ' E~ {x) P ^j , 

where p is the density matrix for the state of the illumi- 
nating field, and is a generalized cross-section for the 
process. This result holds whenever the optical field may 
be considered stationary, quasi-monochromatic and reso- 
nant [HI, Ei| . It can be seen that the absorption process 
will be strongly influenced by the statistical properties of 
the light [1|. 
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B. Master equation and solution using the 
positive-P distribution 



Following a standard model for the NDPO, the Marko- 
vian master equation for the reduced density operator for 
the pump, signal and idler fields is given by [13, Hl[ , 



£_! _ 1. 

dt ih 



Hi,p 



+ 7 



( 2a »/ 5 "l - 4 



aip— pa J Si 



i=i 



+73 (2a 3 p4 - &la 3 p - pala 3 J , 



(2) 



where Hi here excludes the Hi oss term of Eq. dXJ) , which 
couples the cavity modes to reservoir modes. We now 
map this master equation into a classical Fokker- Planck 
equation using the positive-P rep resentation introduced 
by Gardiner and Drummond [15j . For the current prob- 
lem this representation is defined as follows, 



f 3 



(oii*\ a i) 



■P(d) 



(3) 



The six variables oi,oi*, 02,02*, 03 an d a 3* are inde- 
pendent complex variables, and we have written a = 
(oi, Oi*, 02, 02*, 03, 03*). It should be emphasized that 
the asterisks in the variable indices do not correspond 
to complex conjugation. The complex variables on and 
o^ correspond to the mode operators via the relations: 
Q>i\o£i) = oii\oii) and (oi*|a| = o«(o«|. The distribution 
function P(a) may be assumed to have the mathemati- 
cal properties of a probability density function: it is real 
valued, positive, and normalized to one when integrated 
over the full domain D of o. For the master equation, 
Eq. ©, the corresponding positive-P function satisfies 
the following, 



dP 
~dt 



9 , \ / \ 

(701 - K02*03 ) + (702 - KOl*03 1 + 

oai oa 2 

d , s , 9 . . 

(7ai*-K0203*)-|-T; (702*-Kai03*)-|- 

aai* aa 2 * 

9 , 

- — (73^3 + KO1O2 - 7 3 e) + 
oa 3 







da 3 

d 2 a 3 



(7303* + kqi*o 2 * - 7 3 e) + 
3 2 a 3 * 



+ K 



da\da 2 <9ai*<9a2* 



(4) 



This has of the form of a multivariate Fokker-Planck 
equation (HI, and the use of the positive-P distribu- 
tion ensures that the diffusion matrix is positive. The 
Langevin equations corresponding to Eq. Q can be writ- 



ten as, 

da 1 
~df 
da 2 
~dT 
da 3 
dt 

C?Oi* 

dt 
da 2lf 

dt 
da 3 * 

dt 



-701 + K0 2 *0 3 + 



-702 + KOi*0 3 + 



-7303 - KOi0 2 + 7 3 £, 




-701* -I- K0 2 3 * -I- V Ka 3* 
-702* -I- KOl0 3 * + ^/KOs* 



Ki (*)+*& (*)] 
Ki (*)] 

V2 



[& (*)-*& (*)] 
[&(*)+*&(*)] 

V2 



-730-3* 



KOi*0 2 * + 73£, 



(•5) 



where the £i are real-valued, white-noise, stochastic vari- 
ables with (£,-(t)) = and (£i(t)€j(t)) = 5ijS(t - t') for 
all values of the indices. 

Since it is assumed for the mode loss parameters that 
73 3> 7, the pump field can be adiabatically eliminated. 
Setting 0:3 = 03* = then leads to, 

K 

o 3 = e oio 2 , 



03* 



73 

K 

= e 01*02* 

73 



(6) 



Substituting Eq. (J6]) into Eq. 



we obtain, 



da 1 
~dr~ 

da 2 
~~dT 

doi* 
dr 

da 2 * 
~dT~ 
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-OS2- 



-oi*-f 



=— Q2*"f 



(cr-2oiQ 2 ) 2 * / (7 — 1(X\Ol2 ( £l+^2 



n 

(a— 2oio 2 ) oi* 
n 

(a— 2oi*o 2 *) o 2 
n 

(a— 2oi*o 2 *) a.\ 



n 



(7 — 2oi02 



no 



IT— 2oi*02* 



no 



no 



a— 2oi*a2* 
n 



V2 
V2 



(7) 



where time has been scaled in terms of the cavity lifetime 
r = jt. Parameter no = 2773/K 2 is proportional to 
the square of the number of photons in the cavity at 
threshold, and sets the scale for the number of photons 
necessary to explore the nonlinearity of interaction, a = 
273 c/k is a scaled measure of the pump field amplitude 
relative to its value at threshold. At threshold the rates 
7 and eK are equal, and therefore er = n . 

Using the reparametrization set out in Refs. [l2l Il3j|. 
this set of Langevin equations can be decoupled, and the 
dimensionality of the problem can be reduced from eight 
to four. We define four real- valued parameters 771,77,2,773, 
and U4 according to, 
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The Ui may be interpreted as scaled pseudo-quadrature 
variables. An exact distribution for the positive-P func- 
tion for the NDPO may now be written down, as pre- 
sented first in [131, valid for below-, near-, and above- 
threshold regimes. For nn ~ 10 6 — 10 s , a typical value 
for laboratory systems [23], the positive- P function is 
given to a very good approximation by 

P (u) cx exp jd! (ul+ul)+a 2 (1*2+1*4) 



(ul+ul + ul + uf) 2 + 4 (ul + uf) (1*3+1*4) j, 



(9) 



where u— (1*1,1*2,1*3,1*4). Parameters 01, 0,2 correspond 
to the pump strength and are given by 




(10) 

where r = ct/hq = Ke/7. It can be seen from Eq. (|9|) 
that any moment of the form, 



TL\ Tin n3 TlA t-> / -*\ 

u 1 u 2 1*3 U 4 P (it) , 



has the value zero if any of the is odd. The distribu- 
tion of Eq. ([9]) can be used to evaluate the steady-state 
expectation values for any normally-ordered product of 
the field operators. In the next section, we will compute 
the absorption rates for multi-photon recording media, 
exploiting the symmetries of P (u) and making approxi- 
mations valid for the different regimes of operation of the 
NDPO. 



III. 



RESULTS 



In this section, we compute the multi-photon absorp- 
tion rates at the recording medium for a quantum- 
lithographic process, by computing expectation values 
using the positive-P distribution presented in Eq. 
We first propagate the signal and idler variables a,\, ai* 
and «2, «2* through the imaging apparatus, as described 
in Sec. HU By interfering the signal and idler fields at a 
symmetric 50:50 beam splitter, the fields at the output 
are given by 



Pi 
fa 
Pu 
fa* 



{—otx + ia.2) 

(-CK2 + iotl) 

(— au - ia-2*) 

{-0.2* - ia-i*) 
71 ' 



(11) 



Propagating the fields to the multi-photon recording ma- 
terial, the combined field, at a location corresponding 



to an optical phase difference of (f>, is given by a% = 
(fte** + fa) and a 3 * = (ft.e^ + fa*). Substituting 
Eq. (fTT|) in these expressions, we obtain, 

a 3 = J-[ ai (-e^ + i)+a 2 (ie^ -1)] 
v2 

a 3 , = -4 [aui-e-** - i) + a^i-ie-** - I)] .(12) 

The rate of p-photon absorption then is given by the 
average quantity, 



I p 3 {<t>) = (a 3 {uYaUuf 
Substituting Eq. ((8|), we also have, 
1 



(13) 



2c 2 



(1*1-1*2) cos ( - J +i (1*3-1*4) sin ; 



,1*1 + 1*2) cos 



,-1(1*3+1*4) sin } 



(14) 



We first evaluate the fringe pattern for a one-photon ab- 
sorber, given by I 3 ((f)). Inspecting the symmetries of 
the positive-P function, Eq. ([9]), we see that symme- 
tries exist between the variables 1*1 and 1*3, as well as 
1*2 and u 4 , so that P(u\, U2, 1*3, 1*4) = P(l*3, u 2 ,u±, 1x4) 
and P(i*i, u 2 ,U3, 1*4) = P(t*i, 114,1*3,1*2). It follows imme- 
diately that 



Il{4>)=r^{{u\)-{ul)) 



(15) 



and there is no dependence on <f> either below or above 
threshold, and the illumination of the substrate is con- 
stant across the surface in the ensemble-averaged sense. 
Next we consider the case of multi-photon absorbing ma- 
terials. 



A. Below threshold regime 

Inspecting the positive-P distribution, Eq. ©, for the 
below threshold case, nt -C 7, we find that a\ and 0,2 
are large negative quantities, and the distribution can be 
approximated by, 



P(u) & \ —zr 




I 2 

a 2 \u 4 



.(16) 



For a typical NDPO with no ~ 10 6 this approximation is 
very good for the parameter r ranging from to 0.99. The 
four variables Ui are Gaussian and independent. Their 
even order moments are given by 



(uf) = (uf) 
<«?) = «) 



(2k)\ 



k\ (4|ax|) 

(2*0 ! 
k\ (4|a 2 |) A 



k ■ 



(17) 
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and their odd order moments vanish. 

In order to compute the multi-photon absorption rates 
if (0), defined by Eq. (fT3J) and Eq. ([Tg]l. a suitable group- 
ing of terms must be found. To this end we define the 
function, 

F (ui, Uj,4>) = [ui cos (0/2) — iv,j sin (0/2)] 2 , (18) 
so that if can be expressed as 



if(0)=(rV2^) P £ 



X (ULU4, (F (U2,« 3 , </>f~ k ). (19) 

We exploit the fact that the in are not coupled, to evalu- 
ate the moments of separable sums independently. Since 
the pairs of variables u\ and U3, and U2 and U4, share the 
same distribution, we also conclude 



(-ir fe /F(u 2! u3,0r fe 



\p—k 



F(u 3 ,U 2 ,0 + 7r) 

= (F( Ul , U4 ,0 + ^) p - fc ). 
Expanding an arbitrary moment k of F (iti , 1*4, 0) we find 



F(ui,u 4 ,< 



E 

Z=0 



>(fe-0 



x cos 2 ' (0/2) sin 2 ^ 2 ' (0/2) . (20) 



By substituting the expressions for the even moments ui 
and Ui given above, and the form of parameters cii and 
ai defined by Eq. (fTU|) , this sum can be evaluated as, 



^F(ui,U 4 ,0) fc ^): 



ffl [r + cos (0)f . 

fc!4fe(V2^) (l-r 2 ) fe 

(21) 

Using these relations, we arrive at a compact expression 
for the p-photon absorption rate, 




FIG. 2: Plots (a) and (b) show the absorption rates if(0), 
scaled between and 1, for p = 1, • • ■ , 6. It is assumed that 
the NDPO is operated below threshold, with pump param- 
eter r = 0.15 in (a) and r — 0.9 in (b). As p increases the 
corresponding fringe patterns become better defined. 




FIG. 3: The visibility, defined by V = 

[if (</> max )-I$ (0min)] / [ifWmaxH/f (<£min)] and normal- 
ized between and 1, is plotted for p = 2, • • ■ , 6 as a 
function of pump parameter r for the NDPO operating below 
threshold. 



1 p p 



4(l-r 2 ) 



fe=0 



(2fc)! (2p-2fc)! 
/fc! 2 (p-/fc)! 2 



x [r + cos (0)] ' [r — cos (0)] 



(p-fc) 



(22) 



Looking in detail at the form of the absorption rates 
if (0) in Eq. P2^) we see for the following. For the p = 1 
case, the absorption rate is constant and no fringe pat- 
tern is created, as already seen following Eq. (fT5|) . For 
the case of p > 2, if (0) is given as a sum of even pow- 
ers of cos(0). The corresponding fringe pattern therefore 
has terms with periods corresponding to cos(2n</>) where 
2n < p. As p increases, there is a greater contribution 
from higher power terms with sharper interference pat- 
terns. The fringe patterns for the cases of p = 1, • • • , 6 
are plotted in Fig.^a) and (b), for the pump parameter 
r = 0.15 and 0.9 respectively. 



Fig. [3] shows fringe visibilities as r ranges from to 
0.95. It is seen that for every case with p > 2, the vis- 
ibilities are close to 1 when the NDPO is operated far 
below threshold, but as the pump power is increased the 
visibilities fall steadily. In a lithographic process, it is 
possible to compensate for unwanted constant exposure 
by using a substrate with greater depth, and visibilities 
greater than 0.2 are often considered adequate in prac- 
tice. This criterion is satisfied in all cases considered 
here with p > 2. Making a comparison with the case 
of an OPA source reported in the earlier work of Agar- 
wal, et al., in Ref. Q, we find the following. Eq. (|2"2"1) is 
listed explicitly in Table [T] of Appendix [X] for the cases 
of p = I,-- - ,6. Inspecting Eqs. (19-22) of Ref. [1] we 
find the formulae for the multi-photon absorption rates 
identical on identifying r in current analysis and tanh(G) 
in the analysis of Agarwal, et ah, for which G represents 
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the single-pass gain for the OPA. The range for param- 
eter r for the NDPO below threshold (from — > 1_), 
corresponds to values of parameter G for the OPA from 
— * oo. Visibilities for the sub-threshold NDPO are 
listed in Table HT1 of Appendix [Al for p = 1, • • • ,6, which 
agree with the results reported for the OPA across the 
corresponding range for the parameter G. 



component leads to the R(-) contribution, defined by, 
R(S,ai) = f °° p s e-^ p2 - ai ) 2 dp. Writing, 

fl(2«+2t+l,oi) = /' 00 p 2s + 2t - 1 ( P 2 - ai ) e-^-^dp 
Jo 

+ a 1 R(2s + 2t-l,a 1 ), (25) 



B. Near and above threshold regime 

As the NDPO is pumped more strongly and passes 
through threshold, pump rate r exceeds unity, the na- 
ture of the positive- P distribution changes significantly. 
Inspecting the governing parameters a± and a 2 , defined 
by Eq. (|10[) . we see that while a 2 continues to take large 
negatives values of increasing size as r increases, is 
zero at threshold and takes positive values for r > 1. 
The positive-P distribution, given by Eq. ([9]), may now 
be well approximated by the expression below, which is 
valid both near and far above threshold for values of no 
greater than 10 4 (|14j|). 



(23) 




where N denotes a normalization factor for the v.1,113 
component, defined by 



V2 



[7r 3 /2 er fc (- ai /V2)] 



and erfc(z) = (2/y/n) J™ e~ s ds is the complementary 
error function. It can be seen that variables u 2 and U4 are 
independent Gaussian variables with pump parameter a 2 
in the below-threshold regime. However, variables u\ and 
113 with pump parameter ai are now strongly coupled. 
As a consequence, the decomposition given by Eq. (|19p . 
used to evaluate I p (4>) in the below-threshold case, can 
no longer be applied. 

To proceed, in this case we look in detail at the mo- 
ments for the Uj that arise in computing I p (4>) above 
threshold. The moments for the Gaussian variables u 2 
and U4 are as before given by Eq. (|17p . By reparametriz- 
ing m and 113 in polar coordinates, it follows that, 

(ufuf)(ul m )(ul n ) = R(2s + 2t+l, ai ) B (s + l, t + ^j 

(2m)! (2n)\ 



x 2N- 



m\ (4\a 2 \) m n\(4\a 2 \y 



(24) 



B(v) here denotes the Beta function, defined by 
B (S,T) = r(5)T(T) /T(S + T), and arises from the 
integration over the angular component associated with 
u\ and U3. The integration over the remaining radial 



for s + 1 > 1, and integrating by parts, 



p 24+a *- 1 (/j a -ai)e-M^-«i] !, dp 



(if s + t = 1), 



= < 2 v " " ' " (26) 

{(s+t-l)R(2s + 2t-3 1 a 1 ) (ifs + t>l), v ' 

we get a recursion relation for R (•). Here we list first few 
radial functions, 



R ^=2^N 



2 V 2"" 1 J 2nN 
^(5,ai) = ycxp ^-I fll 2 ^ + (1 + ai 2 )^ 

R (7, ax) =(l + ^ exp (-^«i 2 ) + (3oi + a?) ^ 
nt c , ft 1 z\ f 1 a\ (3 + 6 ai 2 + ai 4 ) 



(27) 



Below threshold (r = Ke/7 < I), the parameter a\ is 
negative. It is zero at threshold (r = fte/7 = 1), and 
positive above threshold (r = Ke/7 > 1). As a% increases 
from negative to positive values the NDPO goes through 
a phase transition, and the intensities of the signal and 
idler modes increase very rapidly. For a typical value 
of no = 10 6 — 10 8 , the region for phase transition is very 
narrow and a small change in r from 0.99 — > 1.01 changes 
the statistics of the NDPO dramatically [TEE!- In tllis 
region, we find that the visibility for p-photon absorption 
changes from the below-threshold to the above-threshold 
behavior. 

For the NDPO operating much above threshold 
(ai >> 1), we can considerably simplify the expres- 
sion for If (4>). It follows from Eqs. and J2ZJ) that 
{ufuf) ~ B (s + 1/2, t + 1/2) a\ +t /-K. Since powers of 
no appear only in the denominator for moments of u 2 
and M4, and in the expansion of I p ((j>) the powers of the 
variables u% satisfy s + t + m + n= p , it follows that all 
contributions from variables u 2 and W4 can be neglected. 
Within this approximation, 
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Visibility V 




0.97 0.98 0.99 1 1.01 1.02 1.03 



FIG. 4: The visibility for p = 2 versus pump parameter r as 
operation of the NDPO changes from below threshold to the 
above threshold regime for no = 10 6 . 



and we find the expectation value, 




/ ai rV2^u"\ P ^ (2s)!(2p-2s)! 
V 8 J sl 2 (p-s)\ 2 

x [l + cos(0)] s [l-cos(0)] p ~ s . (28) 

The p-photon absorption rates for much above thresh- 
old, are therefore seen to generate fringe patterns of a 
form independent of the strength of the pump. Fig. [4] 
shows how the fringe visibility changes near the threshold 
regime for the 2-photon case. Note that as r changes from 
0.97 — > 1.03 the parameter eti changes from —42 — » +42. 
As expected, much above threshold the visibility is al- 
most constant. 



IV. DISCUSSION 

In conclusion, we have found that a NDPO source 
can be used to generate fringe patterns with an effective 
wavelength half that for the signal and idler modes which 
interfere at the recording media, when a p-photon absorp- 
tion process is available withp > 2. These fringe patterns 
have high visibility, falling from 1 at low pump power, to 
asymptotic values of 0.2 or greater at high pump power. 
Above threshold, the forms of the fringe patterns are 
insensitive to the strength of the pump, whereas below 
threshold, the forms of the fringe patterns depend on the 
pump power, and tend towards the asymptotic case as 
the NDPO threshold is approached. Comparing with the 
results reported in Ref. [6|], we find that the fringe pat- 
terns generated using an NDPO operating below thresh- 
old (r from — > 1_) are similar to fringe pattern gener- 
ated by an OPA operated with a corresponding gain (G 
from — > oo) . In the latter case there is no cavity and no 



above threshold regime, and the process of optical para- 
metric amplification occurs in propagating modes. This 
result is perhaps surprising, since the earlier analysis of 
Ref. [f| disregards photon losses, and uses a single-mode 
analysis by ignoring all but one down converted spatial 
mode. The NDPO source has some experimental ad- 
vantages. The signal and idler modes at the outputs of 
the NDPO are collimated because of the use of a cav- 
ity, and higher-powers can be obtained than is the case 
for the OPA, both of which are important in light of the 
small cross-sections for typical multi-photon absorption 
processes. Collimated, high-power, outputs also increase 
the speed at which a substrate may be imaged — a criti- 
cal factor in the mass production of, say, computer chips. 
Finally, the reparameterization of the dynamics of the 
NDPO in terms of the four pseudo-quadrature variables 
Ui above, sheds some light on why the change in fringe 
patterns and visibilities are observed. Far-below thresh- 
old, all four of the Ui contribute approximately equally to 
the absorption process, and behave as independent pa- 
rameters. As the pump power is increased, and threshold 
is approached, parameter a± tends to 0. Since ai appears 
in the denominator of each of the moments of variables 
ui and U3, Ui and U3 make a growing contribution com- 
pared to U2 and U4. Above threshold, variables u\ and 113 
make the primary contribution to the absorption process, 
and the effects of U2 and U4 can be neglected. Variables 
u\ and 1*3 are also strongly coupled in this regime. The 
analysis of this paper also contributes compact general 
formulae for the multi-photon absorption rates, Eq. (j2"2")) 
and Eq. fl2SJ). 
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APPENDIX A: 



FRINGE PATTERNS BELOW 
THRESHOLD 



Based on the general formula for the absorption rates 
below threshold, Eq. (|22p . the following tables list explic- 
itly the absorption rates, and the visibilities for the cor- 
responding fringe patterns, for p-photon absorbers with 
p ranging from 1 to 6. 



TABLE I: p-photon absorption rates as a function of the 
optical-path-length difference <j> and the pump parameter r, 
varying from far below to near threshold (from to approxi- 
mately 0.99). 



3r* 



PI 

y [cos 2 (<rt + 2r 2 ] 
[3cos 2 ((£) + 2r 2 ] 



(T^f 

3r z [3 cos 4 (<f>) + 24r 2 cos 2 (0) + 8r 4 ] 



15r 6 



[15 cos 4 (0) + 40r- 2 cos 2 (0) + 8r z 



[5 cos 6 (^>) + 90r 2 cos 4 (<j!>) + 120r 4 cos 2 ((p) + 16r 6 



TABLE II: The visibilities, defined by V = 
[if (</>max)— if(</Wn)] / [J| (0max)+l3 (</>min)] and normal- 
ized between and 1, for the fringe patterns with a 
p— photon absorbing recording material. The third column 
lists the limiting values for the visibilities as the pump 
parameter r ranges from far below to near threshold (from 
to approximately 0.99). 



Visibility V(r) 







1 

1 • 
1 - 



[l+4r 2 ] 

4r 2 
[3 + 4r 2 J 
16r 4 



[3+24r 2 + 16r 4 ] 
1 — 16r 4 

[l5+40r 2 +16r 4 ] 
32r 6 



[5+90r 2 + 120r 4 +32r 6 l 



V(0 



1/(1" 



1 -> 0.20 
1 -» 0.43 
1 -» 0.63 
1 -> 0.77 
1 -> 0.87 
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